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Abstract 

Homogeneous  mixing,  where  all  instances  of  contacts  between 
any  two  members  of  the  population  are  equally  likely,  is  a  common 
assumption  in  modeling  biodefense  policies  against  smallpox.  Such  a 
mixing  pattern  is  rather  unlikely  to  represent  population  interaction  in 
a  modem  urban  setting,  which  typically  is  separated  into  households 
on  one  hand,  and  into  daily  meeting  sites  such  as  schools  and  offices, 
on  the  other  hand.  In  this  paper,  we  develop  a  dynamic  two-level 
social  interaction  model  where  individuals  move  back  and  forth 
between  home  and  daily  meeting  sites,  possibly  passing  through  a 
general  meeting  site  such  as  a  mass  transit  system  or  other  crowded 
areas.  Based  on  this  difference-equations  model,  we  evaluate  the  effect 
of  situational  awareness  and  response  measures,  such  as  vaccination, 
quarantining,  and  movement  restrictions,  on  the  spread  of  smallpox  in 
the  population. 


1.  Introduction 

Responding  to  a  bioterror  attack  of  smallpox  has  become  a  major  concern  to 
governments,  local  public  officials,  and  health  authorities.  This  concern  has  been 
reflected  in  studies  that  model  and  evaluate  possible  response  policies  against 
smallpox  [l]-[8],  A  common  assumption  in  these  models  (e.g.,  [l]-[5])  is 
homogeneous  mixing,  where  all  instances  of  contacts  between  any  two  members  of 
the  population  are  equally  likely.  In  other  words,  interactions  in  the  population  are 
uniformly  random.  The  main  implication  of  this  assumption  is  that  interaction 
probabilities  among  individuals  depend  only  on  the  relative  size  of  the  susceptible  and 
infectious  populations — a  feature  that  simplifies  the  analysis. 

Such  a  mixing  pattern  is  quite  unlikely  to  represent  actual  interactions  in  an  urban 
setting  where  the  population  is  typically  divided  into  interconnecting  subsets. 
Halloran  et  al.  [6]  present  a  heterogeneous  mixing  simulation  model  for  smallpox 
where  ad  hoc  social  structure  is  considered.  The  model  is  applied  to  a  small 
population  of  2,000  people.  An  excellent  review  of  recent  smallpox  models  is  by 
Furgeson  et  al  [7].  A  number  of  studies  examine  nonhomogeneous  mixing  in  other 
epidemic  settings.  Some  ad  hoc  social  mixing  patterns  are  studied  in  [9]  and  [10]. 
Ball  and  Lyne  [11]  consider  a  population  partitioned  into  households,  with  local 
mixing  within  households  and  global  mixing  throughout  the  population,  and  develop  a 
vaccination  optimization  model.  The  effects  of  a  similar  social  structure  are  studied  by 
Koopman  et  al  [12],  Other  heterogeneous  mixing  models  have  been  studied  by  Kaplan 
with  respect  to  the  AIDS  epidemic  [13]  [14]. 

The  concept  of  small  world  networks  [15]  [16]  is  utilized  by  several  researchers  to 
model  nonhomogeneous  transmission  in  a  population  [  1 7]-[  1 9] .  Eubank  et  al.  [18] 
develop  a  detailed  large-scale  urban  traffic  simulation,  and  find  that  interactions 
among  people  form  a  strongly  connected  small-world-like  graph.  They  examine 
several  response  policies  and  conclude  that  outbreaks  can  be  contained  by  a 
combination  of  targeted  vaccination  and  early  detection. 

In  this  paper  we  develop  a  two-level  social  interaction  model  that  consists  of 
households  and  other  daily  meeting  sites  such  as  schools,  offices,  and  mass  transit 
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systems.  This  SIR-based  difference-equations  model  captures  dynamic  features  of 
daily  contacts  among  individuals  in  a  major  urban  area.  We  apply  this  model  to  a 
large  urban  area  (9  million  people)  and  evaluate  the  effect  of  situational  awareness 
(early  detection  and  response)  and  several  response  measures,  such  as  mass 
vaccination,  quarantine,  closure,  mass-transit  shutdown,  and  voluntary  self-quarantine 
on  the  spread  of  the  epidemic  and  on  the  total  number  of  casualties. 

The  rest  of  the  paper  is  organized  as  follows.  Section  2  describes  the  social  structure 
that  forms  the  base  for  our  model  and  analysis.  Section  3  outlines  the  stages  of  the 
epidemic  and  discusses  possible  response  actions.  The  two-level  model  is  described  in 
Section  4.  In  Section  5,  we  report  the  results  of  the  analysis  that  is  based  on  our 
model.  Discussion  and  concluding  remarks  are  presented  in  Section  6.  A  detailed 
description  of  the  difference-equation  model  is  given  in  the  Appendix. 

2.  The  Social  Structure 

We  assume  that  during  each  time  period  (i.e.,  a  day)  a  person  interacts  with  other 
persons  mainly  in  two  places:  in  the  household  (HH)  and  in  the  daily  meeting  site 
(DMS),  such  as  school  or  workplace.  There  may  also  be  incidental  contacts  in  public 
places  such  as  mass  transit  systems,  restaurants,  or  theaters.  We  consider  these 
contacts  as  occurring  in  a  general  meeting  site  (GMS).  During  the  course  of  a  day,  a 
person  is  in  (close)  contact  people  with  a  relatively  small  number  of  individuals  in  the 
HH,  then  she  meets  colleagues,  fellow  students,  or  coworkers  in  the  DMS,  and  finally 
she  may  also  contact  (mostly  strangers)  in  the  GMS.  Figure  1  presents  this 
interaction  pattern. 
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Households 


DMS 


Thus,  we  assume  that  the  population  is  divided  into  m  HHs  of  size  h  each.  There  are 
k  DMSs,  k>  h,  and  one  GMS.  On  each  day,  members  of  a  HH  visit  certain  DMSs.  We 
assume  a  worst-case  scenario,  from  the  epidemic  spread  point  of  view,  where  the 
DMSs  are  chosen  randomly  and  independently  by  HH  members  each  day.  We  also 
assume  that  effectively  no  two  members  of  a  certain  HH  visit  the  same  DMS,  e.g.,  a 
mother  does  not  visit  the  school  (at  least  not  for  a  significant  length  of  time)  and  her 
child  does  not  visit  her  workplace.  Thus,  members  of  the  same  HH  do  not  interact  in  a 
DMS.  They  interact  in  the  HH  and  possibly  in  the  GMS  (e.g.,  if  both  use  the  mass 
transit  system  on  a  certain  day).  Within  each  subset  of  population — HH,  DMS,  and 
GMS — we  assume  homogeneous  mixing.  However,  the  contact  rates  (and  hence  the 
transmission  rates  of  the  disease)  are  different  in  the  three  environments;  the 
transmission  rate  is  highest  in  a  HH  and  lowest  in  the  GMS. 

The  spread  of  the  epidemic  is  observed  at  discrete  time  periods  (days).  Each  time 
period  is  divided  into  two  parts:  the  HH  subperiod  and  the  DMS  subperiod.  During 
the  HH  subperiod  individuals  stay  in  their  respective  HHs  (homes),  while  during  the 
DMS  subperiod  they  are  present  in  their  respective  DMSs  (workplace,  school,  etc.). 
Some  individuals  may  visit  also  the  GMS  in-between  the  two  periods.  The  state  of  the 
epidemic  in  the  HHs  and  DMSs  is  monitored  at  four  points  during  a  time  period:  at 
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the  beginning  and  end  of  the  HH  subperiod,  and  at  the  beginning  and  end  of  the 
DMS  subperiod.  At  the  beginning  of  the  HH  subperiod,  we  observe  the  state  of  the 
members  of  a  HH  after  they  return  from  the  DMSs  and  possibly  GMS,  and  at  the  end 
of  that  period,  we  observe  the  transitions  that  have  occurred  in  the  HH  during  the 
HH  subperiod.  Similar  observations  apply  to  the  DMS  subperiod. 

3.  The  Epidemic  and  Possible  Response  Actions 

The  stages  of  the  epidemic  are  S,  A,  B,  /,  and  Q  for  susceptible,  infected  and 
immunable  (vaccine  sensitive),  infected  and  nonimmunable,  infectious  and 
quarantined.  A  HH  is  said  to  be  infected  if  at  least  one  member  in  the  HH  is  infected, 
but  no  one  is  infectious  (symptomatic).  An  infected  HH  may  be  immunable  (denoted 
type  A)  if  all  of  its  infected  members  are  immunable.  Otherwise,  it  is  not-immunable 
(i type  B).  Clearly,  some  members  in  a  nonimmunable  infected  HH  (type  B)  may  be 
susceptible  (at  stage  S)  or  infected  and  immunable  (stage  A).  A  HH  is  said  to  be 
infectious  ( type  I)  if  at  least  one  member  in  the  HH  is  or  has  been  infectious,  and  it  is 
said  to  be  quarantined  ( type  Q )  if  it  has  been  put  in  quarantine.  Otherwise,  a  HH  is 
said  to  be  susceptible  ( type  S).  We  assume  that  vaccination  and  quarantine  are  applied 
to  HHs  and  not  to  individuals.  We  assume  perfect  vaccination  efficacy,  therefore, 
susceptible  or  infected-immunable  HHs  that  are  vaccinated  are  removed  from  further 
consideration.  Once  an  infectious  individual  is  detected,  his  entire  HH  is  quarantined. 
If  that  HH  has  not  been  previously  vaccinated,  all  asymptomatic  members  are 
vaccinated  upon  entering  the  quarantine.  Only  infectious  HHs  are  quarantined. 

HHs  of  types  B  and  /  that  are  vaccinated  are  labeled  BV  and  IV,  respectively.  Also, 
QV  denotes  a  quarantined  HH  that  has  been  previously  vaccinated.  Figure  2  presents 
the  transitions  among  the  stages. 


4 


Figure  2:  The  Stages  of  the  Epidemic 

Since  we  assume  perfect  vaccination  efficacy,  vaccinated  HHs  of  type  S  and  A  are 
immune  and  removed.  Households  of  types  B  and  /  may  be  vaccinated,  but  only  the 
individuals  at  stages  S  and  A  in  those  HHs  become  immune.  The  rest  (those  at  stages 
B  and  I)  are  unaffected  by  the  vaccination. 

Without  loss  of  generality  we  assume  that  transitions  in  the  stage  of  a  HH,  including 
vaccination  and  quarantining,  occur  during  the  HH  period.  A  DMS  may  be  either 
infectious,  if  at  least  one  member  is  infectious,  or  noninfectious.  Infectious  DMSs  may 
generate  new  infected  individuals. 

We  consider  the  following  response  actions: 

•  Mass  vaccination; 

•  Quarantine; 

•  Shutdown  of  GMSs  (e.g.,  shutdown  of  a  mass  transit  system); 

•  Closure  of  DMSs  (e.g.,  closing  up  schools);  and 

•  Encouraging  people  to  stay  home. 
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The  response  actions  are  initiated  after  a  certain  number  of  individuals  become 
infectious  (Stage  /).  We  assume  that  the  vaccination  and  the  DMS  closure  processes 
start  after  there  are  A1  infectious  individuals.  This  number  is  the  vaccination/DMS 

closure  threshold.  The  smaller  the  value  of  Aj ,  the  better  the  situational  awareness 

capabilities  of  the  system.  The  shutdown  of  the  GMS  is  also  triggered  by  the  number 
of  infectious  individuals  in  the  population.  The  GMS  is  shut  down  if  this  number 
exceeds  A, .  While  the  GMS  is  shut  down  instantaneously  following  a  decision  to  that 
effect,  the  process  of  closing  up  the  DMSs  is  gradual  and  takes  time.  Once  the  DMS 
closure  process  is  initiated,  it  proceeds  with  a  rate  8 . 

We  assume  that  a  fraction  y  of  the  population  passes  through  the  GMS,  and  a 
proportion  f3  of  the  (not  yet  isolated)  infectious  individuals  complies  with  requests  of 

the  authorities  and  voluntarily  stay  at  home.  Note  that  in  terms  of  the  epidemic  spread, 
the  latter  assumption  represents  a  worst-case  scenario;  self-imposed  quarantine  does 
not  apply  to  susceptible  individuals — only  to  infectious  ones.  Finally,  infectious  HHs 
are  isolated  (removed)  at  a  rate  p . 

4.  The  Model 

The  model  comprises  a  set  of  difference-equations  shown  in  the  Appendix.  It  is 
essentially  a  deterministic  model,  but  with  embedded  probabilities  that  capture  the 
transition  intensities  between  stages.  The  spread  of  the  disease  is  observed  at 
two  levels:  high  level,  at  which  we  observe  the  transitions  between  sets  (HHs  and 
DMSs),  and  low  level,  at  which  we  observe  the  transitions  within  sets.  At  any  time 
period  t,  we  record  the  number  of  sets  of  a  certain  type  (e.g.,  S,  A,  etc.)  at  the  high 
level,  and  the  average  profile  (composition)  of  disease  stages  within  a  set,  at  the  low 
level.  For  example,  B(fi)  is  the  number  of  HHs  of  type  B,  and  sB(t),aB(t),bB(t )  are 

the  average  numbers  of  individuals  in  such  HHs  that  are  at  stages  S,  A,  and  B, 
respectively.  The  difference  equations,  shown  in  the  Appendix,  describe  the 
transitions  between  sets  of  stages  at  the  high  level,  and  the  changes  in  the  average 
profiles  within  sets,  at  the  low  level.  Recall  that  at  each  time  period  (day)  the 
epidemic  is  observed  four  times:  at  the  beginning  and  end  of  the  HH  subperiod,  and 
the  beginning  and  end  of  the  DMS  subperiod. 
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To  simplify  notation,  we  use  capitals  to  denote  both  the  type  of  a  set  and  the  number 
of  sets  of  that  type.  For  example,  I(t )  denotes  an  infectious  HH  and  the  number  of 
such  HHs  in  the  population.  This  double  meaning  should  cause  no  confusion. 

Specifically,  the  symbol  XJ  ( t )  denotes  the  number  of  sets  (HHs  or  DMSs)  of  type  X 
at  time  t.  The  index  j  is  0, 1  where  j  =  0  indicates  a  beginning  of  a  subperiod  (HH  or 
DMS),  and  j  =  1  indicates  the  end  of  it.  To  summarize,  X  gets  the  following  values: 

S  -  Number  of  susceptible  HHs. 

A  -  Number  of  infective-immunable  HHs. 

B  -  Number  of  infective  nonimmunable  HHs. 

BV -  Number  of  infective  nonimmunable  HHs  that  have  been  vaccinated 
(Only  individuals  at  stage  B  remain  infective,  the  rest — S  and  A  individuals — 
are  vaccinated  and  removed). 

/  -  Number  of  infectious  HHs  that  have  not  been  vaccinated  yet.  fy  are  newly 
infected  HHs. 

IV-  Number  of  infectious  HHs  that  have  been  vaccinated.  VIq  are  newly 
infected  HHs. 

Q  -  Number  of  isolated  HHs. 

QV -  Number  of  isolated,  previously  vaccinated,  HHs. 

D  -  Number  of  open  DMSs. 

ID  -  Number  of  open  infectious  DMSs. 

The  notation  at  the  low  level  is  of  the  form  yJx  (/) ,  where  y  indicates  the  stage 

epidemic,  X  is  the  type  of  HH  or  DMS,  and  j  is  a  0,1  parameter  as  before, 
for  example: 

Average  number  of  susceptible  individuals,  at  the  beginning  of  the  /-th  HH 
period,  in  an  infective-immunable  HH  that  has  not  been  vaccinated  yet. 
s\v(t)-  Average  number  of  susceptible  individuals,  at  the  end  of  the  t- th  HH  period, 
in  an  infectious  HH  that  has  been  vaccinated. 
a°A(t)~  Average  number  of  immunable  infective  individuals,  at  the  beginning  of  the 
t- th  HH  period,  in  an  infective  HH  that  has  not  been  vaccinated  yet. 


of  the 
Thus, 
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b)n(t)-  Average  number  of  nonimmunable  infective  individuals,  at  the  end  of  the 
t- th  HH  period,  in  an  infective  nonimmunable  HH  that  has  been  vaccinated. 
i\{t)~  Average  number  of  infectious  individuals,  at  the  end  of  the  t-  th  HH  period, 
in  an  infectious  HH  that  has  not  been  vaccinated  yet. 

S°fD(t) ~  Average  number  of  susceptible  individuals  in  an  infectious  DMS  at  the 
beginning  of  the  DMS  period. 

In  addition,  we  denote 

a^w(t)  -  Average  number  of  newly  infected  individuals  in  an  infectious  DMS,  who 

belong  to  a  HH  of  type  X,  X=S,  A,  B,  BV,  /,  IV  (t)  is  the  average 
total  number  of  newly  infected  individuals  in  an  infectious  DMS. 

Note  that  while  the  number  of  individuals  in  a  HH  remains  constant  throughout  the 
epidemic,  the  average  number  of  individuals  in  a  DMS  changes  over  time  as  HHs  are 
isolated  and  infectious  persons  stay  home. 

The  parameters  of  the  model  are: 

M  -  Number  of  HHs. 

K  -  Number  of  DMSs. 
h  -  Size  of  a  HH. 

V  -  Vaccination  capacity  (#  of  HHs/day). 

V  -  Vaccination  rate  (percentage  of  population  vaccinated). 
p  -  Transition  rate  between  stages  A  and  B. 

q  -  Transition  rate  between  stages  B  and  I. 
aH  -  Infection  rate  in  a  HH. 

CCD  -  Infection  rate  in  a  DMS. 

CCG  -  Infection  rate  in  the  GMS. 

y  -  Fraction  of  individuals  that  visit  the  GMS  during  a  time-period. 
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P  -  Fraction  of  infectious  people  in  a  (not  yet  isolated)  infectious  HH  that  stay 
home  because  of  personal  choice  (regardless  if  their  corresponding  DMSs  are 
closed  or  open). 

p  -  Isolation  rate  of  infectious  HHs. 

G  -  Recovery  rate  from  isolation. 

S  -  Closure  rate  of  (infectious)  DMSs. 

Aj  -  Vaccination/DMS  closure  threshold. 

A2  -  GMS  shutdown  threshold. 

To  demonstrate  the  basic  idea  of  the  model  we  present  next  a  sample  of  three 
typical  equations: 

1 .  High  level  transition  during  the  HH  subperiod: 

Bl(t)  =  [B°  (0(1  -  qfe(t)  +  A0  (0(1  -  (1  -  p)a° 0)  )]Max{  1  -  v(t),  0} .  ( 1 ) 


Equation  (1)  gives  the  number  of  HHs  of  type  B  (not  yet  vaccinated)  at  the  end  of  the 
HH  subperiod. 

2.  Low  level  transition: 


4(0  =  —  slG)B\m-qf«Hs\(t)A\t){\-{\-Py>)  Max{l-v(0,0}  .  (2) 


B\t)l 


Equation  (2)  gives  the  average  number  of  susceptible  individuals  in  HHs  of  type  B. 

3.  Low  level  in  a  DMS. 


(3) 

Equation  (3)  gives  the  average  number  of  newly  infected  in  an  infectious  DMS. 

The  complete  set  of  equations  is  presented  in  the  Appendix. 


New/.\  0  -0 

aID  V  )  -  aDsIDlID  • 
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5.  Analysis 

We  consider  two  base  cases.  In  Base  Case  1,  we  assume  that  the  response  policy  is 
based  only  on  mass  vaccination  and  quarantine  of  infectious  HHs.  Other  response 
measures  such  as  DMSs  closure,  GMS  shutdown,  and  compliance  with  self-imposed 
quarantine  are  not  implemented.  Base  case  2  is  the  complement  of  Base  Case  1.  There 
is  no  mass  vaccination,  only  quarantine  (and  vaccination)  of  infectious  (type  I)  HHs 
that  are  detected.  However,  DMSs  are  gradually  closed  up,  the  GMS  is  shut  down 
after  a  while,  and  a  certain  proportion  of  infectious  individuals  (not  yet  quarantined) 
stay  home.  In  both  cases,  we  assume  that  the  initial  attack  resulted  in  ten  casualties 
(infected  people)  in  each  one  of  five  DMSs  plus  ten  infected  in  the  GMS.  Table  1 
presents  the  parameters  that  are  fixed  for  both  base  cases  and  in  the  subsequent 
sensitivity  analysis. 


Parameter 

Description 

Value 

M 

Number  of  HHs 

3,000,000 

h 

Average  size  of  a  HH 

3 

k 

Number  of  DMSs 

10000 

aH 

Infection  rate  in  a  HH 

0.5 

°D 

Infection  rate  in  a  DMS 

0.001 

aG 

Infection  rate  in  the  GMS 

0.000001 

P 

Disease  stage  A  rate 

0.3 

q 

Disease  stage  B  rate 

0.12 

p 

Disease  stage  /  rate 
(Isolation  rate  of  infectious  HHs) 

0.3 

0 

Disease  stage  Q  rate 
(Recovery  rate) 

0.083 

r 

Proportion  of  population  that  passes  through  the  GMS 

0.5 

Table  1:  Values  of  Fixed  Parameters 


Table  2  presents  the  values  of  the  remaining  parameters  in  the  two  base  cases. 


Parameter 

Description 

Base  Case  1 

Base  Case  2 

5 

Closure  rate  of  DMSs 

0 

.03 

P 

Fraction  of  infectious  individuals  that  stay  home 

0 

.5 

A, 

Vaccination/DMS  Closure  Threshold 

20 

20 

a2 

GMS  Shutdown  Threshold 

No  Closure 

250 

V 

Vaccination  Capacity  (HHs/Day) 

157,000 

0 

Table  2:  Values  of  Variable  Parameters 
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Based  on  these  parameters,  Base  Case  1  (vaccination,  no  public  restrictions)  results  in 
2,068  infectious  individuals  in  addition  to  the  casualties  of  the  initial  attack,  while 
Base  Case  2  (no  vaccination,  public  restrictions)  results  in  2,070  additional  casualties. 
Evidently,  the  total  numbers  of  casualties  are  essentially  equal.  That  is,  preventive 
measures  that  include  closure  of  DMSs  at  a  rate  of  3%  per  day,  shutting  down  the 
GMS  when  there  are  250  infectious  cases,  and  50%  infectious  “stay-home” 
compliance  is  equivalent  to  the  mass  vaccination  of  157,000  HHs  per  day  with  no 
other  social  movement-control  measures.  However,  the  way  the  epidemic  evolves 
over  time  in  these  two  cases  is  significantly  different,  as  shown  in  Figure  3. 


—  Base  Case  1 
Base  Case  2 


Time  (Days) 


Figure  3:  Daily  Numbers  of  Newly  Infectious 

Figure  4  presents  the  daily  number  of  people  in  quarantine  in  both  base  cases. 


Time  (Days) 


—  Base  Case  1 
Base  Case  2 


Figure  4:  Daily  Numbers  of  People  in  Quarantine 
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As  shown  in  Figure  3,  the  epidemic  in  Base  Case  1  is  shorter,  but  with  a  higher  daily 
peak  than  Base  Case  2.  Figure  4  shows  the  ramification  of  this  effect,  namely,  higher 
demand  in  Base  Case  1  for  peak  quarantine  capacity  than  in  Base  Case  2. 

Next,  we  investigate  the  sensitivity  of  these  base  cases  to  changes  in  the  values  of  the 
operational  parameters.  The  measure  of  effectiveness  is  the  number  of  casualties, 
which  are  the  total  number  of  infectious  individuals.  Figure  5  shows  the  impact  of 
situational  awareness  and  responsiveness  of  the  public-health  system.  The  delay  in  the 
vaccination  process,  which  is  represented  by  the  number  of  infectious  individuals  that 
trigger  the  initiation  of  the  process,  affects  the  number  of  casualties.  If  the 
public-health  system  has  excellent  situational  awareness  and  is  quick  to  respond,  then 
the  minimum  possible  number  of  casualties  is  about  850.  If,  because  of  poor 
situational  awareness  and/or  slow  response,  there  are,  say,  80  infectious  people  in  the 
population  before  vaccination  starts,  then  the  number  of  casualties  is  higher  in  an 
order  of  magnitude. 

Figures  6-8  depict  the  consequences  of  Base  Case  1  (with  V  =  150,000  HHs/day) 
when  additional  social  movement-control  measures  are  imposed.  Figure  6  investigates 
the  effect  of  DMS  closure  on  the  number  of  casualties.  With  no  DMS  closure,  the 
number  of  casualties  is  close  to  2,400.  A  closure  process  at  a  rate  of  10%  a  day 
decreases  the  number  of  casualties  to  about  500.  The  effect  of  the  delay  in  shutting 
down  the  GMS  is  shown  in  Figure  7.  If,  theoretically,  the  GMS  (e.g.,  mass  transit 
system)  is  shut  down  immediately  following  the  attack,  then  there  are  only 
282  casualties.  If  it  is  shut  down  when  there  are  already  500  infectious  individuals, 
then  the  consequence  is  1,987  casualties. 
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Figure  5:  The  Effect  of  Vaccination  Threshold  -  Base  Case  1 


DMSs'  Closure  Rate 


Figure  6:  The  Effect  of  DMSs  Closure  Rate  -  Base  Case  1 


Figure  7:  The  Effect  of  GMS  Shutdown  Threshold  -  Base  Case  1 
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The  effect  of  voluntary  self-quarantine  is  demonstrated  in  Figure  8.  Note  that  if  60% 
of  the  infectious  individuals  stay  home,  then  the  number  of  casualties  is  reduced  by 
almost  80%  compared  to  no  self-quarantine. 


Figure  8:  The  Effect  of  Self-Quarantine  -  Base  Case  1 

Let  us  consider  now  Base  Case  2  and  investigate  the  effect  of  adding  to  the  social 
movement-control  measures  mass  vaccination.  Figures  9-12  depict  the  effect  of 
vaccination  capacity,  with  respect  to  DMSs’  closure  rate  ( 8 ),  vaccination/DMS 
closure  threshold  ( A, ),  GMS  shutdown  threshold  ( A2 ),  and  self-quarantine 
compliance  rate  ( /? ),  respectively.  Each  chart  represents  the  base  case  and  two 
additional  possible  scenarios,  worse  case  (darker  shade)  and  better  case 
(lighter  shade). 


Daily  Vaccination  Capacity 


- 5  =.01 

- 5  =  03 

5  =  05 


Figure  9:  Effect  of  DMSs’  Closure  Rate 
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Daily  Vaccination  Capacity 


Figure  10:  Effect  of  Vaccination/DMS  Closure  Threshold 


- A2=500 

- A2=250 

A2=50 


Daily  Vaccination  Capacity 


Figure  11:  Effect  of  GMS  Shutdown  Threshold 
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Daily  Vaccination  Capacity 


Figure  12:  Effect  of  Self-quarantine 


Figures  9-12  reveal  a  consistent  picture,  namely,  that  the  impact  of  situational 
awareness  and  social  movement-control  measures  can  be  significant,  but  it  decreases 
as  the  vaccination  capacity  increases.  One  small  exception  is  the  vaccination/closure 
threshold  Aj .  Even  if  the  vaccination  capacity  is  240,000  HHs/day,  the  number  of 

casualties  in  the  worse  case  scenario  A,  =50  is  almost  3.5  times  higher  than  in  the 
case  where  A1  =5. 

6.  Discussion  and  Conclusions 

Figures  3  and  4  demonstrate  that  (mass)  vaccination-free  policies,  where  the  response 
relies  entirely  on  social  movement-control  measures,  may  be  compatible  to  a  policy 
that  is  based  on  mass  vaccination.  Moreover,  while  the  latter  policy  eradicates  the 
epidemic  faster  than  the  former,  the  required  peak  quarantine  capacity  is  significantly 
higher.  Figures  5-8  show  what  happens  if  a  mass  vaccination  policy  is  augmented 
with  social  movement-control  measures  and  situational  awareness  is  considered.  The 
effects  are  quite  significant.  For  example,  if  one-third  of  the  infectious  population 
withdraws  to  their  home,  the  number  of  casualties  is  reduced  by  55%  compared  with 
no  self-quarantine.  A  similar  effect  is  recorded  when  the  GMS  is  shutdown  early  in 
the  epidemic.  If  the  trigger  for  a  shutdown  is  250  infectious  people  in  the  population, 
then  the  number  of  casualties  is  1,594.  If  the  trigger  is  50  people,  then  the  result  is 
only  742  casualties. 

These  results  indicate  that  a  hybrid  policy  that  combines  moderate  rate  of 
mass-vaccination  with  moderate  application  of  social  movement-control  measures 
may  be  an  efficient  response  policy.  Table  3  presents  the  parameters  of  a  third 
base-case,  which  is  a  combination  of  the  two  approaches:  vaccination  and  social 
movement-control. 


Parameter 

Description 

Base  Case  3 

5 

Closure  rate  of  DMSs 

.01 

P 

Fraction  of  infectious  individuals  that  stay  home 

.3 

A, 

Vaccination/DMS  Closure  Threshold 

10 

a2 

GMS  Shutdown  Threshold 

100 

V 

Vaccination  Capacity  (HHs/Day) 

75,000 

Table  3:  Base  Case  3 
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The  total  number  of  casualties  in  Base  Case  3  is  1,358 — about  35%  less  casualties 
than  in  Base  Cases  1  and  2.  Similar  to  Figure  3,  Figure  13  compares  the  three  base 
cases  over  the  period  of  the  epidemic. 


Figure  13:  Daily  Numbers  of  Newly  Infectious  -  Base  Cases  1,  2,  and  3 


While  the  peak  of  the  epidemic  is  about  the  same  as  in  Base  Case  2  (47  casualties  in 
Base  Case  2,  43  casualties  in  Base  Case  3),  the  epidemic  is  eradicated  much  faster: 
76  days  compared  with  1 13  days. 
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Appendix:  Difference-Equations  Model 


A.  HH  Subperiod 

Let 


v(0  = 


m 

S°(t)  +  A°(t)  +  I°(t) 


(4) 


v(t)  is  the  vaccination  rate. 


High  Level 


S\t)  =  S°(t)Max{\-v(t),Qi} 

(5) 

A1  ( t )  =  A0  (0(1  -  p)a°  {t)  Max{  1  -  v(0, 0} 

(6) 

Bl  (t)  =  [5°  (0(1  -  qfe(t)  +  A0  (0(1  -  (1  -  p)a°(0  )]Max{l  -  v(t),  0} 

(7) 

BVX(t) 

=  BV°(t)(l-q)b*At)  +  [B°(t)(\-q)bAA +A°(t)(\-(\-p)aAA)]Min{v(t),\}  (8) 

4  (0  =  B°  (0[1  -  (1  -  qfAt)]Max{  1  -  v(t),  0} 

(9) 

IVq  ( t )  =  5F°[1  -  (1  -  q)b-  {t)]  +  B° (0[1  -  (1  -  q  fM)]Min{\,  v(t)} 

(10) 

l\t)  =  1°  (t)Max{  1  -  v(0, 0}  (1  -  p) 

(11) 

IV 1  (0  =  (IV°  (t)  + 1°  (t)Min{\,  v(0})(l  -  P) 

(12) 

Q(t)  =  I°(t)p  +  Q(t- 1) 

(13) 

QV(t)  =  IV°(t)p  +  QV(t- 1). 

(14) 

Low  Level 

Susceptible  HH  ( S) 

s1s(t)  =  s°s(t)  =  h 

(15) 

a\(t )  =  b\(t)  =  4(0  =  a°s(t)  =  b°s(t)  =  4(0  =  0 . 

(16) 
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Infective  Immunable  HH  (A) 


s\  (t)  =  s°A(t) 

(17) 

aA  (0  =  «S(0 

(18) 

b]4(t)  =  b°A(t)  =  0 

(19) 

^  — 

II 

^  O 

II 

O 

(20) 

Vaccinated  Infective  Not-immunable  HH  (B  V) 

4r  (0  =  S°BV  (0  =  aBV  (0  =  «8|/  (0  =  llv  (0  =  4f  (0  =  0  (25) 

*;,(o = — )tt  + 

o/7  (0L 

“(^(0  +  a»5°(0(i-^(o+  .  (26) 

(0 - —  (0(1  -  (1  -  pfA (n)  Min{v(t),  1} 

1-(1  -p)aAt) 
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Newly  Infectious  HH  ( Iq ) 


(0=4(0 

(27) 

a)(t)  =  a)B{t)(\-  p) 

(28) 

f  \ 

bl»  (t)  =  bB(t)  1  1  +aB(t)p 

{  1-(1  -qf-m) 

(29) 

i  (0  =  62(0 - 

(30) 

Vaccinated  Newly  Infectious  HH  ( IVo ) 

siv0  (0  —  siv0  (0  ~  aiv0  (0  —  fl/F0  (0  —  6 

(31) 

iff  ) 

h\rU)’-  .  />2,-«)  1 - SF0(<)(1-(1-«)V(0) 

N(0  \  {  1  - (1  -qf"' 10  J 

+  *J(0  1 - ^7-  5°(0(l-(l-^(°)Mm{v(0,l}  (32) 

l  1-(1  -q)K{t))  \ 

40(0  =  -\—\qb*v(t)BVQ(t)  +  qb  °B(t)B°(t)Min{v(t),l}} 

IVq  ft)  ■ 

Infectious  HH  (I) 

s)  =  sl(^~aH^l) 
a\ft )  =  cij(t)(l-  p)  +  aHSjij 
b)(t)  =  b(f(t)(\-q)  +  a(j(t)p 
ilj(t)  =  ij(t)  +  bj(t)q . 

Vaccinated  Infectious  HH  ( IV) 


s\yff)  ~  siv  (0  ~  a/v(0  -  aiv( 0  _  6 

(38) 

^/F  (0  =  biv  (0(1  ~d) 

(39) 

oj 

+ 

II 

(40) 
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(34) 

(35) 

(36) 

(37) 


Isolated  Not  Previously  Vaccinated  HH  (Q) 

(Assumption:  Individuals  not  previously  vaccinated  are  vaccinated  immediately  upon 


arrival  at  the  quarantine). 

SQ(t)  =  SQ(t)  =  ag(t)  =  ag(t)  =  0  (41) 

6e (<> = {('■ -  NqO -  mt- - 1}  ■ +  (!  ■ -  « )  bf  (0/°  (.Dp)  (42) 

iQ{t)  =  iQ{t-\)  +  ^[qbg{t-\)Q{t-\)  +  {qbQI{t)  +  i0im\t)p  }.  (43) 

Isolated  Previously  Vaccinated  HH  ( (Q  V) 

sqv(  0  =  sQy(t)  =  aQv(0  =  aQv({) =  0  (44) 

V(0  -  -  DO  -  dQVO  -  0  +  *  m-Nv°(»p}  (45) 


iQV{t)  =  iQvit-\)  +  -P—{qbQV(t-l)QV{t-\)  +  (qb,j¥(t)  +  i,>Ir(t)W\t)p}.  (46) 
At  the  end  of  the  HH-cycle,  the  total  number  of  susceptibles  and  immunable 


infected  are: 

4otai  (0 = s1  m + a1  (t)s\  (t)+B[  (t)slB  0 1 ) + 4  (o40  (0 + (04  (0  (47) 

40fa/  (0 = 41  (o4  (o + b{  (04  (0 + 4  <7)4,  (o + 71  (o4  (0  •  (48) 

And  the  total  number  of  infectious  individuals  is: 

Lai  (0 = 4  o40  (0) + 71  (0h  (0) + (o40  (0) + ivl  (o4  (0)  •  (49) 

The  number  of  commuting  infectious  individuals  from  an  infectious  HH  is 

ic\  (t)  =  (1  -  f3)i)  (t)  .  (50) 

ic\  (t),ic}y(t)  and  ic,v^  (t)  are  defined  similarly. 

The  total  commuting  infectious  individuals  is: 

Lotalit)  =  4  (0*4.  (0)  +  71  {t)ic\  (0)  +  IVi  (t)ic)Vo  (0)  +  IV1  (t)icjv  (t)  .  (51) 
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B.  Transition  HH  -»  GMS  ->DMS 


DMS 


D(t)  =  (l-S)D(t-l) 


(52) 


ID(t)  =  D(t) 


1- 


1- 


*  y 


l- 


fcj(ov'(0' 

K 


1- 


■  i 

lciv\t) 


K 


J 


•  l 


1  — 


zc 


7K, 


(0 


N^oCO^ 


a: 


(53) 


Explanation:  Suppose  D(t)  =  uK.  Only  a  proportion  u  of  the  DMSs  are  open  and  therefore  only  a 
proportion  u  of  the  population  leaves  home.  Since  each  member  of  a  HH  goes  to  a  different  DMS,  this 
means  that  only  ui  leave  home.  ui/uK  =  i/K. 

Individuals 


We  assume  that  the  contacts  in  the  GMS  occur  between  the  HH  cycle  and  the 
DMS  cycle. 


'CMS  (0=7  — ^  icTotal  0  )  (34) 

sGMS  (0  =  Y^Totali^  (35) 

The  number  of  newly  infected  individuals  at  the  GMS  is: 

aGMS  (0  =  aGsGMSkjMS  •  (56) 


Let 

New 


a, 


GS 


(0.47(0. 


47(0and47(0, 


47(0 + 47(0 + 47(0 + 47(0 = 4S  (0. 


denote  the  number  of  newly  infected  at  the  GMS  that  belong  to  S,  A,  B,  and  /  HH, 
respectively.  Since  the  probability  that  a  newly  infective  belongs  to  a  certain  type  of  a 
HH  is  proportional  to  the  number  of  susceptibles  in  such  a  HH,  we  have: 


New s,\  New  /,*.  S  (f)h 

aGS  \l>-aGMS\l)  i 

sTotal(t) 

(57) 

„New „New  (t)SA(t) 

aGA  \l>-aGMSyl>  1 

sTolal(f) 

(58) 

„NeWf+\  „New  (O^gO) 

"C5  \l>-aGMS\l>  i  .  . 

sTotal(t) 

(59) 

New s.\  New  ^  ^SI0  (0  +  ^  0)5/  (0 

«G/  (0  =  aGMs(t)  ]  _ 

5 Total  (0 

(60) 
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Infectious  DMS  (ID) 


,o  .  _  D(t)  icTolal(t) 

lID  0)~ 


K  ID(t ) 


,0  sf  \  _  sTotal  (0  “  «GM5  (? ) 
sIDO  )  ~  ~ 


(61) 

(62) 


Since  response  actions  and  transitions  between  disease  stages  are  assumed  to  take 
place  during  the  HHs  cycle,  we  do  not  need  to  track  either  the  noninfectious  open 
DMSs  or  the  individuals  who  are  at  the  latent  stages  (immunable  and  nonimmunable) 
of  the  disease. 


C.  DMS  Cycle 

Since  transitions  between  disease  stages  are  assumed  to  take  place  only  in  the  HHs, 
the  DMSs  do  not  change  their  status  during  this  cycle.  The  only  parameter  of  interest 
during  the  DMS  cycle  is  the  number  of  newly  infected. 

Individuals 


Infectious  Open  DMS  (ID) 

The  number  of  newly  infected  individuals  in  an  open  DMS  is: 

News.\  0  -0 

aID  (t )  “  aDsIDlID  ■  (63) 


A  newly  infected  individual  may  belong  to  a  susceptible  HH  (S),  an  infective  HH  (A), 
or  an  infectious  HH  (I). 


Let 

New 


a 


DS 


s.\  New/-.\  News.\ _ j  News.\  Newr.\  ,  News.\  ,  New,.\  ,  News.\  News.\ 

0),aDA  (0; aDB  (t)  ar|d  aDI  (t),  aDS  (t)  +  aDA  (t)  +  aDB  (t)  +  aDS  (t)-aID  (t), 


denote  the  number  of  newly  infectives  that  belong  to  S,  A,  and  I  HH,  respectively. 
Since  the  probability  that  a  newly  infective  belongs  to  a  certain  type  of  HH  is 
proportional  to  the  number  of  susceptibles  in  such  a  HH,  we  have: 


New/.\  News.\  $  f)h 

aDS  >  ~  aID  {U  i 

sTotal(t ) 

(64) 

„New „A rewf+\B 
aDB  \l)  ~  aID  V)  i 
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D.  DMS  -  HH  Transition 
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(t)  =  A1  (t  - 1)  +  7  (t  - 1)  -  S°  (0 
=  1) 


(68) 

(69) 

(70) 
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(71) 
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Susceptible  HH  ( S ) 
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Infective  Not-immunable  HH  (B) 


^(o=^[51a-i)4a-i)-^-iKr(^-i)-«Gra-1)]  (so) 

a°B(t)  =  [B1  it  - 1)4 {t - 1)  +  m  - 1  )aND7 (t- 1)  +  aNJ (?-  !)]  (81) 

b°B(t)  =  bl(t- 1)  (82) 

i°B(t)  =  0.  (83) 


Infectious  HH  (I) 
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